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Abstract 

Experiments aiming at measuring the neutron electric dipole moment (nEDM) 
are at the forefront of precision measurements and demand instrumentation of 
increasing sensitivity and reliability. In this paper, we report on the develop- 
ment of a dedicated acquisition and control electronics board for the nEDM 
experiment at the Paul Schcrrer Institute (PSI) in Switzerland. This multifunc- 
tion module is based on a FPGA (Field-programmable gate array) which allows 
an optimal combination of versatility and evolution capacities. 



1. Introduction 

The search for a non-zero electric dipole moment (EDM) of the neutron, 
initiated 60 years ago by Norman Ramsey, continues today to motivate exper- 
imental activity with increasing precision. A non-zero electric dipole moment 
for a spin 1/2 particle, such as the neutron, implies the violation of both par- 
ity (P) and time reversal (T) symmetries; according to the CPT theorem this 
would also imply a violation of the combined CP symmetry. Since CP violation 
beyond the standard model of particle physics is required to explain the asym- 
metry between matter and antimatter, the so-called baryon asymmetry of our 
Universe, electric dipole moments provide a window on the early Universe when 
the baryon asymmetry was generated. For a recent review on the connection 
between neutron physics and cosmology, see [1] and references therein. 

The best measurement of the neutron EDM [2] was obtained at the Institute 
Laue-Langevin (ILL) in Grenoble, France and provided the upper limit |d n | < 
2.9 x 10~ 26 e.cm (90% C.L.). It was extracted from the precession frequency of 
trapped ultra-cold neutrons (UCN) in a weak (Bo= 1 uT) magnetic field, using 
Ramsey's method of separated oscillating hclds. An atomic magnetometer using 
a vapor of polarized 199 Hg atoms cohabiting in the same volume as the UCN was 
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used to control and correct for the magnetic field fluctuations. The experiment, 
installed at the Paul Scherrer Institute (PSI), Villigen, Switzerland, uses an 
upgraded version of this apparatus and aims at improving the sensitivity by an 
order of magnitude. 

The development of a completely renewed data-acquisition system (DAQ) 
was part of the numerous upgrades. To this end, new electronic modules were 
needed to execute all requested actions and control the experiment. Instead 
of multiple specialized modules, we have opted for a central and multifunction 
module. With such a solution, we expect to benefit from the limited number of 
connectors, thus minimizing the delicate problem of bad contacts. 

After a brief recall of the experiment principle, a description of a typical 
measurement cycle will be given in section 2. Section 3 will present in details 
the electronic board and section 4, the embedded firmware. Finally, section 5 
will report on some of the achieved performances. 

2. Description of the measurement sequence 

The PSI EDM experiment uses the RAL/SUSSEX/ILL apparatus, con- 
nected to the newly built PSI UCN source. The principal elements of the 
apparatus are depicted in fig. 1. The central part, the precession chamber, 
is a 21 liters cylindrical bottle used to store polarized ultra-cold neutrons and 
polarized 199 Hg atoms. This chamber is shielded against the ambient magnetic 
field with a four layer mu-metal assembly. A highly homogeneous vertical mag- 
netic field Bo= 1 pT is generated inside the mu-metal shield by a set of coils. 
Additionally, oscillating transverse (to Bo) fields can be applied by pairs of coils 
in the X or Y directions (fig. 1). The bottom and top parts of the precession 
chamber also serve as electrodes to generate a strong electric field (typically 
120kV/12cm), either parallel or antiparallcl to the Bo magnetic field. 

The principle of the experiment consists in measuring the neutrons Larmor 
precession frequency using Ramsey's method of separated oscillating fields. The 
neutron Larmor frequency is proportional to the magnitude of Bo , plus possibly 
a small frequency shift due to the neutron EDM. A frequency difference, pro- 
portional to E, between the parallel and antiparallel configuration would be the 
signature for a non-zero EDM. However, the small drifts of the Bo field during 
the experiment could easily hide the signal. To measure and correct for these 
fluctuations, a mercury comagnetometer is used [8]: the precession frequency of 
mercury atoms, also proportional to the Bo magnitude, occupying the same vol- 
ume as UCN - hence the name comagnetometer - is measured via the intensity 
modulation of a resonant light beam traversing the storage chamber. 

A 'standard EDM cycle', sketched in fig. 2, is composed of seven successive 
steps detailed below. 

1. UCN fill: Polarized UCN are guided toward the storage chamber. The 
filling time, controlled by the UCN valve, lasts typically 20 s. 

2. Mercury fill: To let the 199 Hg gas, optically pumped in the polarisation 
cell, go in the precession chamber, the mercury valve is opened for 2 s. 
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Figure 1: Experiment description. 



3. Mercury pulse: A rotating transverse magnetic field is generated for 
2 s by two pairs of coils. The signal frequency is equal to the Larmor 
frequency of the "Hg atoms (7.59 Hz for Bo=luT), measured during 
the previous cycle. The signals feeding the X and Y coils are in phase 
quadrature and their amplitudes are adjusted so that the Hg spins rotate 
by 7r/2rad. A triangular envelope was chosen in order to minimize the 
effect of the mercury pulse on the neutrons spins. 

4. First UCN pulse: As for the mercury pulse, two pairs of coils are used 
to generate for 2s a rotating field inducing a tt/2 flip of the neutrons 
spins. The signal frequency is also calculated from the mercury Larmor 
frequency measured at the previous cycle. However, for the neutrons, we 
use four working points whose frequencies arc close to the neutron Larmor 
frequency (29.16Hz for B =lpT). 

5. Free precession: For a typical duration of 200 s, both neutron and mer- 
cury spins precess in the horizontal plane around Bo- The mercury free 
precession is optically monitored, using a through-going polarized light 
produced by a mercury lamp. The light intensity, measured with a photo- 
multiplier (PMT), is modulated at the Larmor frequency of the mercury 
spins. The PMT signal is sampled with the DAQ electronics and an es- 
timate of the frequency is extracted. As discussed above, the estimated 
frequency /n g serves to set the frequency of both mercury and neutron 
pulses at the next cycle, thus accounting for the drifts in magnetic field. 

6. Second UCN pulse: At the end of the precession time, a second neutron 
7r/2 pulse, in phase with the first one, is generated, thus completing the 
Ramsey sequence. 
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Figure 2: Snapshot of 'standard EDM cycle'. 



7. Neutron polarization measurement: The UCN valve is opened to 
empty the precession chamber and neutrons fall down toward a neutron 
detector. On their way, neutrons encounter a magnetized ferromagnetic 
foil that serves as spin analyzer by letting only one spin component go 
through. Once all neutrons with the right spin component have been 
counted, an Adiabatic Fast Passage (AFP) spin flipper is turned on to 
count the number of neutrons in the other component. From this infor- 
mation and the knowledge of the n/2 pulse frequency, the neutron Larmor 
frequency can be extracted. 

The described sequence is repeated continuously to accumulate statistics, 
with periodic reversal of the electric field. Two years of data-taking at the PSI 
UCN source are planned to record about 50,000 cycles and gain an order of 
magnitude in sensitivity. 

3. Hardware description 

The electronic- board block diagram is shown in fig. 3. It is designed around a 
Field Programmable Gate Array (FPGA), whose main purposes are to perform 
the system precise operation sequencing ('micro-timer') and the data acquisi- 
tion. The communication with the board is done via a USB2 capable micro- 
controller, this allows to perform the system control and the data readout. 
Furthermore, it is used to load the FPGA firmware at board power up and to 
perform the monitoring of 16 input channels with a 16 bit ADC. It directly con- 
trols the signal multiplexer and reads out the ADC. The channel is reserved 
for reading the Photo Multiplier Tube (PMT) DC signal. 

A Rb atomic clock, operating at 10 MHz, is directly supplied to the FPGA in 
order to reference the whole system (signal generation and measurements). The 
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Figure 3: Block diagram of the electronic board. 



various sine wave generation is ensured by two kinds of Direct Digital Synthesis 
(DDS) circuits whose outputs are connected or disconnected by FPGA driven 
switches. The first DDS kind (AD9833) is used to provide the spin flipper signals 
(typical frequency of 19 kHz). They have a 28 bit resolution phase accumulator 
and thus can achieve a frequency resolution of 37mHz when supplied with a 
10 MHz reference clock. These DDS are configured by the micro-controller via 
a Serial Peripheral Interface (SPI) link. Each DDS output signal, which has an 
amplitude resolution of 10 bit, is filtered by an active second order Low Pass 
Filter (LPF) before passing through the FPGA controlled switch. The switch 
output is then amplified by 16 before being provided to the spin flipper coils. 

The second DDS kind (AD9852) is used to provide the two neutron (around 
30 Hz) and two Mercury (around 8 Hz) excitation signals with a very high pre- 
cision in frequency. They feature a 48 bit phase accumulator and given the 
fact that they are referenced at 500 kHz (a divided version of the atomic clock), 
they have a tuning resolution of 1.77 nHz. Also these DDS offer the possibility 
to add an offset to the phase accumulators result and given the fact that each 
DDS pair is updated synchronously, it is possible to generate the excitation 
signal pair with a known phase relationship. Moreover, by directly controlling 
the DDS parameter update by the FPGA, before closing the output switch, it 
is possible to adjust the start phases of both sine waves. Another advantage, 
is that the sine wave output amplitude, which has an amplitude resolution of 
12 bit, can be digitally adjusted either by a constant factor or by using the 
Output Shape Keying (OSK) feature. As illustrated in fig. 4, this allows the 



control of the signal amplitude as a function of time and therefore to generate 
the required triangular shape envelope instead of the abrupt rectangular shape. 
With the DDS clock running at 500 kHz, the ramping time can be adjusted from 
32.768 ms up to 2.096 s with a step resolution of 8.192 ms. The ramping (up or 
down) is detcrministically activated by the micro-timer embedded in the FPGA 
via a dedicated pin (labeled OSK). An example of the control sequence allowing 
to adjust the initial phase, the ramping and the switch control is displayed in 
fig. 5. Each DDS output signal is amplified with an adjustable gain up to 
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Figure 4: Sketch of the output amplitude as a function of time for the abrupt and OSK modes. 
While the abrupt mode toggle the output amplitude instantly to full scale, the OSK feature 
allow the control the sine/cosine amplitude as a function of the time. 
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Figure 5: Example of a full control sequence allowing to adjust the start phase, ramping and 
switch control with a deterministic timing. 



a factor ten and filtered by an active fourth order Low Pass Filter (LPF) be- 
fore passing through the FPGA controlled switch. A serial resistor (820O) has 
been added in the circuit output in order to have a current source like output. 
Consequently, the magnetic field generated is less sensitive to the coil resistance 
variation due to thermal fluctuations (typical coil resistance is around 10 O). 

The PMT signal is connected to a conditioning block detailed in fig. 6. The 
first stage converts the current signal in a voltage signal with a current pream- 
plifier having a 1 Mf2 resistor feedback. The alternating (AC) and continuous 
(DC) components of the signal are then separated: the DC component is ex- 
tracted via a second order low pass filter while the AC part passes through a 
Q = 8 band pass filter centered around 8 Hz. The resulting signal is level shifted 
by half of the ADC full-scale in order to fit in the available ADC range. 

Some additional features of this module include: 

• the possibility to control eight external electromechanical devices syn- 
chronously to the micro-timer with the provided buffered TTL outputs. 
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Figure 6: Block diagram of the PMT signal conditioning. This block converts the current 
signal in a voltage signal and separate the DC and AC parts of the signal. 



• the generation of a real-time four bit signal at each micro-timer step, en- 
coding the step number to permit the synchronization with other electronic 
modules. 

• the availability of twelve 32 bit scalers to record the neutron counts at the 
end of a EDM cycle. 



4. Firmware description 
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Figure 7: Overview of the firmware. It is composed of seven different blocks: the frequency 
divider, the USB interface for communication, the scalers, the micro-timer for precise operation 
sequencing, the ADC interface for data acquisition and the two DDS interface blocks. 



An overview of the FPGA firmware which is the core of the system and whose 
main purposes are to perform the system precise operation sequencing and the 
data acquisition is shown in fig. 7. It is composed of seven different blocks: the 
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frequency divider, the USB interface for communication, the scalers, the micro- 
timer for precise operation sequencing the ADC interface for data acquisition 
and the two DDS interface blocks. 

The frequency divider is used to divide the atomic clock frequency by a 
factor 20 for referencing of the DDS. As explained in section 3, the frequency 
must be lowered in order to reach the highest resolution on the output signal 
frequency but also to be able to have slow ramping amplitudes. This division is 
performed by a to 9 counter which keeps a 50% duty cycle. 

The micro-timer block is designed to manage up to 16 time sequences. Each 
sequence has a programmable duration, that can be adjusted with a 1 us step 
resolution in a 1 to 2 32 -l range, and can manage up to 32 different actions. 
The actions to be executed during each step are controlled by a 32-bit word, 
each action having a dedicated bit. The 32 micro-timer parameters, the 16 
durations and the 16 action masks, are written in the dual port block RAM 
by the USB interface before starting the time sequence execution. The micro- 
timer block is managed by a simple two states Finite State Machine (FSM) 
which starts the time sequence execution when the run signal, provided by 
the USB interface, is activated. Before executing each instruction cycle, the 
corresponding parameters are read from the memory. Additionally, the 4 Most 
Significant Bit (MSB) of the memory read pointer are used to reflect the current 
sequence in progress. This information is available for USB readout and on the 
board for external synchronization. The possible micro-timer actions are listed 
below: 

• The eight 'TTLout' may activate external electro-mechanical actuators 
(UCN and Hg valves). 

• 'HgExcUpd' and 'NeuExcUpd': signals to trigger the DDS parameters 
update for mercury and neutron, respectively. 

• 'enHgExc' and 'enNeuExc': enabling of signals generation for mercury and 
neutron respectively. 

• 'OskHg' and 'OskNeu': initialization of the ramping of the sine wave am- 
plitude (sec figure 4). 

• 'cnAdc': enabling of PMT signal acquisition. 

• The four 'enSpinFlipper' are used to activate the selected 19 kHz sine wave 
generator. 

• The 'cnSpinUpCnt' and 'cnSpinDnCnt' bit are respectively used to allow 
the 12 associated scalers to count during spin up and spin down. 

The DDS interface blocks are used to simultaneously configure the DDS 
pairs (mercury or neutron) with the parameters provided by the USB interface 
in their embedded block memory. Each DDS requires up to 40 configuration 
bytes. Consequently, the two DDS parameters bytes are concatenated in a single 
word and the whole configuration is conveniently stored in the 40 first memory 
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words (16 bit). The FSM, which manages the communication protocol with 
the DDS, initiates the parameter transfer upon reception of the update signal 
provided by the micro-timer table or by the USB interface. 

The firmware features a scaler containing two banks of twelve 32-bit counters, 
each counter having the capacity to count with a frequency up to 100 MHz. 
These counters are operated in a gated mode with the gating signals provided 
by the micro-timer as previously described. 

The PMT signal acquisition is managed by the ADC interface block which 
performs the data acquisition at the selected sampling period and when ac- 
tivated by the micro-timer. The digitized data are written in the embedded 
output FIFO and real time data readout is realized via the USB interface. The 
sampling period, which is a sub-division of the atomic reference clock, can be 
adjusted by steps of 100 ns between a minimum of 10 us and a maximum of 
1677 ms. It should be noted that running at a conversion rate significantly 
higher than the signal frequency (8 Hz) can be an asset if one wants to apply 
an off-line sharp digital filtering to further reduce the bandwidth given by the 
analog band pass filter. 

5. Board performances 

A series of tests was conducted to check the main performances of the board 
and ensure it fulfills the requirements. As an example of the various basic 
functional tests, we show in fig. 8 the generation of a sine wave signal with a 
triangular envelope. This 8 Hz sine wave having a duration of 2 s and 1 s rise 
and fall times, was sampled at 1kHz with the board's own ADC. 

The most critical and demanding function of the board concerns the PMT 
signal treatment. Indeed, as discussed in section 2, the Hg co-magnetometer is a 
central piece of the nEDM apparatus since, by measuring the mercury frequency, 
we can precisely measure the magnetic field seen by the neutrons and correct 
for its fluctuations. Systematic tests were therefore carried out to study the 
measurement chain performances. 

At first, the AC input band pass filter was characterized. To this end, one 
of the mercury '8 Hz' outputs was connected to the PMT signal input, and its 
frequency was varied from 0.2 Hz to 11.5 Hz. The result of the measurement, 
displayed in fig. 9, is in accordance with a Q = 8 filter, centered at 8 Hz. 

The next crucial test consisted in checking that the intrinsic board noise 
was negligible as compared to the physical noise which is obviously one of the 
parameter limiting the precision. The Power Spectrum Density (PSD) of the 
complete board input chain was measured: this includes the amplifier (tested 
with a gain 1), the band pass filter and the ADC. As in the previous measure- 
ment, the input signal was a mercury '8 Hz' signal generated by the board itself. 
A sampling rate of 100 Hz and an acquisition time of 100 s were chosen as they 
represent typical values used by the experiment. The Fast Fourier Transform of 
the digitized signal is shown in fig. 10. The measured noise floor is compatible 
with the expected ADC performances as advertised by the manufacturer for the 
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Figure 8: Capture of an 8 Hz excitation sine wave, having a duration of 2 s and 1 s rise and 
fall time. The signal was sampled at 1kHz with the boards own ADC. 



same amount of points. It is also worth noting that the noise floor inside and 
outside of the filter bandwidth is identical, indicating that the electronics noise 
is due to quantization only. The signal to noise ratio obtained from the PSD 
spectrum is about 14000. 

Those performances should be compared with the precession signal measure- 
ment in the experimental setup. The PSD of a typical measurement is shown 
in fig. 11. In contrast to the loop- back measurement, the noise floor greatly 
increases in the filter bandwidth, showing that the measurement noise is dom- 
inated by the physical noise. The signal to noise ratio obtained from the PSD 
spectrum is about 840. 

From this test, we can conclude that the sensitivity in the frequency extrac- 
tion of the mercury signal will not be limited by the electronic performances 
(security coefficient of 10). 



6. Summary 

Wc have designed and constructed an electronic board (fig. 12) for the data 
acquisition and control of the neutron electric dipole moment experiment at the 
Paul Scherrer Institute. The board, organized around a Field Programmable 
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Figure 9: PMT input bandpass filter frequency response measurement. The filter parameters 
are verified, i.e center frequency lying at 8 Hz and a quality factor of 8. 



Gate Array (FPGA), is a multifunction module which fulfills a large fraction of 
the requirements regarding the time-sequencing and the data acquisition of the 
experiment. Thanks to the use of a FPGA, it is compact, versatile and evolutive. 
Additionally, the high integration level and the corresponding limited number 
of connections should result in a better reliability. A list of the main functions 
and specifications is given in table 1. 
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Figure 12: nEDM experiment acquisition and control electronic board picture. The board has 
a dimension of 300 mm X 150 mm. 
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Function 


specifications 


microtimer 


16 time sequences with 32 actions 
1 us time resolution 


Excitation pulse 


frequency resolution 1.77 nHz 
phase resolution 0.011 ° 
amplitude resolution 12 bit 
output current range ±20 mA 

ramping time: 32.768ms to 2.096s (8.192ms resolution) 


Spin nipper 


frequency resolution 37mHz 
amplitude resolution 10 bit 
output voltage range ±5 V 


PMT acquisition 


16 bit digitization for 3.3 V 

input current conversion factor: 1 V/uA 

bandpass input filter 4 th order, f c =8, Q=8 

selectable voltage gain: 1 to 1000 

sampling rate: 10 ps to 1677 ms (100 ns resolution) 

SNR = 14,000 at full scale 


Slow monitoring 


16 bit digitization 

Input voltage range: 0-3.3 V 

sampling rate of up to 50 Hz 


Neutron counter 


maximun counting rate of 150 MHz 
full scale 2 32 



Tabic 1: Recapitulation of the module specifications. 
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